Fermionic description of spin-gap states of antiferromagnetic Heisenberg ladders 

in a magnetic field 
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Spin gaps— the energy gaps in magnetic excitation 
spectra— of low-dimensional Heisenberg antiferromagnets 
have been attracting considerable interest in recent years. 
This fascinating subject was initiated by Haldane re- 
newed at the opportunity of high-temperature supercon- 
ductivity being discovered ||, and further developed via 
the synthesis of ladder materials ||Q. Dagotto, Riera 
and Scalapino || pioneeringly pointed out a possible elec- 
tronic mechanism for the formation of spin gaps in two- 
leg ladders. Following investigations (6) revealed that as- 
sembling chains into ladders, the crossover between one 
and two dimensions is far from smooth— Heisenberg lad- 
ders with an even number of legs have a spin gap, while 
those with an odd number of legs have no gaps. 

Spin ladders in a magnetic field provide further inter- 
esting topics. Their ground-state magnetization curves 
have extensively been studied |7|-|l0f in an attempt to 
find quantized plateaux. Spin-S* L-\eg Heisenberg ladders 

N / L L-1 \ 
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may exhibit magnetization plateaux at M/M sat = m/LS 
provided LS - in 6 Z [0JT|,|l|], where M and M sat are 
the magnetization and its saturated value, respectively. 
Since coupled-chain materials are likely to be systemati- 
cally obtained in the case of spin i Jl3| , most of numerical 
efforts are being devoted to S = | ladders. Then, with- 
out any bond polymerization, the most tractable system 
of our interest is the three-leg ladder. Lanczos diago- 
nalization of finite systems 0, a series-expansion tech- 
nique H and density-matrix renormalization-group cal- 
culations Jio[ | are all in agreement to support the exis- 
tence of the plateau at M/M sat = 1/3 for strong inter- 
chain coupling J' > J. As for more-than-three-leg lad- 
ders, the existence of plateaux as well as their surviving 
region is still left to be verified. 

Spin-i two- and three-leg ladders are indeed realized 
in layer materials Sr„_iCu n +i02n (n = 3, 5, •••) P,^3|, 
while four- and five-leg ones in La4+4 n Cu8+2nOi4+8n 
(n = 2,3, ••■) 0. It is unfortunate that nevertheless 
large degrees of freedom prevent us from making numeri- 
cal access to multi-leg ladders. In such circumstances, we 
propose a systematic approach to Heisenberg ladders in 
order to describe their spin-gap states in a field. Employ- 
ing the Jordan- Wigner transformation on a unique path 



and then making a mean-field treatment of the fermionic 
Hamiltonian, we semiquantitatively visualize the appear- 
ance of plateaux and estimate the corresponding critical 
interactions J' c at an arbitrary number of legs, L. 

It is along a snake-like path, (n, I) = (1, 1) — > (1, 2) — > 

> (1,L) -> (2,L) - (2,L-1) -f •■• -> (2,1) 

(3, 1) — + • • ■, that we define spinless fermions. This 
elaborately ordered path was first proposed by Dai 
and Su |lq] and turned out to describe the spin gap 
as a function of L much better than the naively or- 
dered path fig] , which is usually employed. When 
we introduce renumbered spin operators S n j = S n> i 
(Sn,L-i+i) for an odd (even) n, the spinless fermions 
are created as c^, = S+iexp[-m(J2"~i J2j=iS^S^ + 

Yli=\ ^n.i^n. ■))]■ Making a standard mean-field treat- 
ment [l5| l of the fermionic Hamiltonian and assuming the 
spatial homogeneity (S* 
we obtain 
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where 



= A^ 1/2 £„e- ifcn 4, 



The effective 



Hamiltonian (H), together with the Zeeman term 

— H *}2 k Y^t=ii^k i°k,i ~ 1/2), is numerically diagonal- 
ized, adopting the open boundary condition along rungs, 
while taking the thermodynamic limit under the periodic 
boundary condition along legs. 

In order to verify the reliability of the present ap- 
proach, we plot in Fig. [l] the thus-obtained spin gaps 
(H c i) and saturation fields (H c2 ), the former of which 
are compared with highly precise numerical estimates 
jTrfl , whereas the latter of which with the exact solu- 
tions H c2 = 2 J + J'[l + cos(n/L)]. Although the present 
scheme somewhat overestimates H c \ with increasing L, 
it correctly tells whether the gap survives or not. As 



1 



for H C 2, the present calculation can be regarded as ex- 
act. We further show in Fig. || typical calculations of 
the ground-state magnetization. All the plateaux satis- 
fying the criterion L/2 — M/N € Z indeed appear with 
increasing J'. Mean- field approaches generally under- 
estimate quantum fluctuations and therefore necessarily 
overestimate the magnetization. Allow us, however, to 
stress that the length of a plateau is still well describable 
in our treatment, which is essential to estimate the lower 
boundaries of the plateau-surviving region. 
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FIG. 1. Spin gaps H c \ (a) and saturation fields H C 2 (b) for 
the I/-leg ladders with J' = J, where the present calculations 
(o) are compared with density-matrix renormalization-group 
estimates [17] (x) in (a) and with the exact values (x) in (b). 



five-leg ladder material (La2Cu04)3La2Cu4C>7 0] en- 
courages us to make theoretical explorations into the 
unique way from one- to two-dimensional antiferromag- 
nets in a field. It is also important, further numerical 
verification of the plateau at M/M sat = 3/5 more sur- 
viving than that at M/M sat = 1/5 for L = 5. 




FIG. 3. Critical interchain interaction to induce magneti- 
zation plateaux, J' c , for the L-leg ladders, where the corre- 
sponding magnetization values M/M sa t are indicated beside 
symbols. 
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FIG. 2. Ground-state magnetization curves of the three- 
and five-leg ladders with J' = 3 J (a) and J' = 5 J (b), respec- 
tively. Numerical-diagonalization results [7] are also shown by 
a dotted line in (a). 

Now we explore the main issue in Fig. || The 
critical value J' c = 1.04 ± 0.01 for L = 3 is in good 
agreement with the previous estimate J' c ~ J |g[ ob- 
tained through the series expansion from the strong- 
rung-coupling limit J/ J' — > 0. J' c appears to increase 
with L, where even- and odd-L ladders may form dis- 
tinguishable series of their own. Here is a conclusive 
remark: Among possible nontrivial quantized magneti- 
zations, M/M sat = 1 - 2/L, 1 - 4/L, • • • ,2/L (1/L), the 
inner plateaux are easier to induce. The plateaux at the 
end-value magnetizations M/M sat = 1 — 2/L,2/L (1/L) 
are generally induced with larger J' than those at the 
inner-value magnetizations are. Practical observation of 
multi-plateau magnetization curves may less be feasible 
with increasing L. At the least, however, an available 
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